Theory of Andreev reflection in junctions with iron-based High-T c superconductors 
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We construct a theory for low-energy quantum transport in normal | superconductor junctions involving the re- 
cently discovered iron-based high-T c superconductors. We properly take into account both Andreev bound 
surface states and the complex Fermi surface topology in our approach, and investigate the signatures of the 
possible order parameter symmetries for the FeAs-lattice. Our results could be helpful in determining the sym- 
metry of the superconducting state in the iron-pnicitide superconductors. 
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Introduction. Very recently, a family of iron-based super- 
conductors with high transition temperatures was discovered, 
with a concomitant avalanche of both experimental and theo- 
retical activity , 'i 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 1 1 » 1 2 . The highest T c measured so 
far in this class of materials is 55 K, and many experimental 
reports indicate signatures of unconventional superconduct- 
ing pairing. However, it remains to be clarified what the exact 
symmetry is for both the orbital- and spin-part of the Cooper 
pair wavefunction - there has for instance been reports of both 
nodal 11 and fully gapped 10 order parameters (OPs) in the lit- 
erature up to now. 

Probing the low-energy quantum transport properties of 
superconducting materials has proven itself as a highly 
useful tool to access information about the symmetry 
of the superconducting OP. 13 The conductance spectra of 
normal | superconductor (N|S) junctions often contains impor- 
tant and clear signatures of the orbital structure of the OP. For 
instance, when the OP contains nodes in the tunneling direc- 
tion with a sign-change across the nodes on each side of the 
Fermi surface, the conductance will display a large zero-bias 
conductance peak (ZBCP) due to the presence of Andreev sur- 
face bound states J 4 - 

Two recent studies 1011 utilized the method of point-contact 
spectroscopy in order to study the symmetry of the supercon- 
ducting OP in LaOo.gFo.i-sFeAs and SmFeAsO0.85F0.15> re- 
spectively. The findings were in stark contrast. Namely, the 
large ZBCP found in LaOrj.gFn.i-aFeAs gave evidence of a 
nodal order parameter, while the data of SmFeAsO0.85F0.15 
clearly indicated a nodeless OP. In both of these studies, 
the Blonder- Tinkham-Klapwijk 15 (BTK) framework was used 
to analyze the data theoretically, using the extension to 
anisotropic pairing by Tanaka and Kashiwaya 16 . In this 
model, one considers a cyndrical or spherical Fermi surface 
with a free-electron dispersion relation, which does not ac- 
count for the non-trivial multiband Fermi surface topology 
and dispersion relation in the iron-pnicitides. One might argue 
that the extended BTK model nevertheless may suffice to de- 
scribe the transport properties of these materials qualitatively, 
but this statement clearly warrants a detailed investigation. 

In this Rapid Communication, we construct a theory of 
low-energy quantum transport properties of the iron-based 
high-Tc superconductors by considering a N|S junction rel- 
evant for point-contact spectroscopy and scanning-tunneling- 
microscopy measurements. In doing so, we model fairly accu- 



rately the Fermi surface topology and the associated quasipar- 
ticle dispersions, in order to see how this affects the results as 
compared with the usual BTK-paradigm. We consider several 
possible OP symmetries which may be realized in the iron- 
pnicitides. 

Theory and results. We adopt the minimal two-band model 
derived in Ref. 8 (Fig. QJ) , in which the normal-state Hamil- 
tonian reads 



(1) 



where the fermion basis 4>^ a — [dk xcr , dk ya ] T contains the an- 
nihilation operators for electrons in the d xz - and e^-orbitals 
with spin a and wavevector k, respectively. We have also de- 
fined £kx = —2t\Cx - 2t 2 Cy - 4:t 3 C X Cy, C^y = ~U^S X S yi 

£kt/ = ~2t2C x — 2t\c y — 4:t3C x c y , with Cj = cos(fcja), Sj = 
s'm(kja), j — x, y, and a is the lattice constant. By diagonal- 
izing the above Hamiltonian, one obtains 



fif = (ekx + e ky )/2 -n± ^(e* x - e ky ) 2 /4 + e 2 xy (2) 

where the new basis ^ a — [7jti>7k ( i] T consists of new 
fermion quasiparticle operators in the bands + and — which 
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FIG. 1: (Color online) (a) Illustration of the two-dimensional FeAs- 
plane with the d xz - and d yz -orbitals and hopping between them, 
as proposed in Ref. 8 . (b) Sketch of the Fermi surface topology for 
the long-lived quasiparticle excitations in a minimal two-band model 
(see main text for parameter values). 
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are hybrids of the d xz - and c?j, 2 -orbitals. The Fermi surface 
topology is given by fl^ = 0, and gives an electron-like band 
(+) and hole-like band (— ) shown in Fig. [TJ? for the choice 
h = -1, < 2 = 1.3, * 3 = U = -0.85, fi = 1.54, all mea- 
sured in units of \t\ |. Our choice of parameter set is motivated 
by the fact that it reproduces the same Fermi surface structure 
as LDA band structure calculations 17 , and was also employed 
in Ref, 1 ^ 7 . The new fermion operators are related to the old 
basis fa a by 



Ck = ekatf/[(eka - £ky)/2 + {e kx - e ky ) 2 / 4 + 4z y L 

i -Ck 



9 k CT ^k 



- 1\ 



, p k = (i + CkV 1/2 x 



Ck l 



(3) 



We now introduce a superconducting pairing between the 
long-lived quasiparticles 7^, A = ±, which then automat- 
ically accounts for both inter- and intra-band pairing in the 
original fermion basis fa: 



Hsc = 



kA 



kp 



h.c. 



(4) 



In this way, we may diagonalize the total Hamiltonian H = 
Hm + Hsc by introducing a final fermion basis rj^ = 



[Ckt; c -kJ T describing the quasiparticle excitations in the su- 
perconducting state. After discarding unimportant constants, 



we find that H = Ek.A^k^L) 1 ^ 



-'ken 



Ei 



A\2 
kJ 



|A^| 2 ] 172 . This result is formally identical to a two-band su- 
perconductor with gaps A£ and normal-state dispersions f^, 
A = ±. The belonging wavefunctions which describe the 
quasiparticle excitations read 

= {[u^e-^fe 1 ^, [ifr^.tiipe-**-'}, 



K) 2 = i - K) 2 - 5 (i + V^ 2 - l A kl 2 /^), (5) 

for quasiparticles with positive excitation energies E > 0. 
Here, k A denotes the Fermi momentum for band A while 



|A£| 



We have now effectively described the superconducting 
state as a two-band model with gaps A k and normal-state dis- 
persions fij. This has allowed us to obtain a simple form for 
the wavefunctions in Eq. (0 that are to be used in the scatter- 
ing problem below. The trade-off for this advantage, however, 
is that the k-dependence of the gap functions A k in general 
will become quite complicated. To see this, we may transform 
Eq. (0| back to the original fermion basis fa by means of our 
expression for Pk in Eq. ^ to find that: 



ky| 



+ Akx^LZ-kyi - di y1 d\ xl ) + h.c.j , (6) 

where Ak x and Ak a are the intra-orbital gaps while Ak xy is 
the inter-orbital gap, defined as A^ x — (A k + Ck^kV^k > 



t J 



A ky = (A k + CkA+)/i/+ , A k „, = Ck(A+ - A k with 
v£ = (1 ± (k). We see that the inter-orbital pairing vanishes 



in the case where A k 



A^. However, we emphasise that 
our model does account for inter-orbital pairing Ak xy , an d 
that Akxy whenever A k ^ A k . Assuming spin-singlet 
and even-frequency pairing, there are three possible s-wave 
symmetries {Aq, Aq(c x + c y ), Aqc x c v } and two possible d- 



wave symmetries {Aq(c 



/), Aos x s y } for the supercon- 



ducting order parameters A^ x and Ak y in terms of the square 
lattice harmonics. The gaps in the ± quasiparticle hybridized 
bands are then obtained as A k = (Akz — CkAkj/)/^ an( ^ 
A k = (Aky — CkAkx)/ 1 ^- Note that the extended s- 
wave symmetry ~ c x c y changes sign on the electron- and 
hole-Fermi surfaces, similarly to the s±-scenario suggested 
in Ref A 

We are now in a position to evaluate the conductance of the 
system. The presence of a Fermi-vector mismatch between 
the normal and superconducting side of the junction is as- 
sumed to be manifested through an effective decrease in the 
junction transmission. Since the Fermi velocity may be dif- 
ferent in the two bands with normal-state dispersions fijr, we 
allow for different barrier parameters Z in the two bands. 
For a specified pairing symmetry, there are then four fitting 
parameters present: the barrier strength Z x and gap magni- 
tude A% for band A = ±. By generalizing the results of 
Refs. ll5lll6l to a two-band model which also takes into account 
the non-trivial Fermi surface topology in Fig. [ij, we obtain 
the following expression for the normalized tunneling con- 
ductance: G(eV)/G = Ex, ky f(k y )a x (eV)/[2f(k y )a x N }, 
where = [1 + (Z A ) 2 ]^ 1 and 



a x (eV) = [a x N [l + a x N \T x + (k, eV)\ 2 + (a x N - 1) 

x|r A (k, e y)r A (k, e y)| 2 ]]/f|i + K - 

xT x + (k,eV)T x Jk,eV)p x (k)\ 2 



1) 



I±(k,eV) 



P x (k) = 



eV- ^{eVy -\A(±\k Xl k y )\* 

\A x (±\k x ,k y )\ 
A x (-\k x ,k v )[A x (\k x ,k y )}* 
\A x (-Xk x ,k y )A x (Xk X7 k y )\ ' 



(7) 



where f(k y ) = cos(k y a/2) is a weighting function that mod- 
els the directional dependence of the incoming quasiparticles. 
The strategy is now to sum the conductance over allowed val- 
ues k y G [— n/a, n/a] for the electron- (A = 1) and hole-like 
(A = —1) Fermi surfaces, and solve for k x from Eq. (fj) by 
fi k = for a given k y . In what follows, we choose equal value 
for the barrier tranparencies Z + = Z~ = Z and gap magni- 



tudes Aq 



A = Aq in the two bands for simplicity, and 



add a small imaginary number 5 to the quasiparticle energy to 
model inelastic scattering: eV — > eV + iS, S/Aq = 1CP 2 . 
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FIG. 2: (Color online) Plot of the conductance spectra for tunneling along the (lOO)-axis in an iron-pnicitide N|S junction for several possible 
order parameter symmetries. Only in the d xy -wscve case A = Aos x s y is there a considerable ZBCP. Note that different scale on the voltage- 
axis for this case due to the narrowness of the ZBCP. High (low) values of the parameter Z denotes low (high) transmissivity interfaces. 



As in Ref. 18 , we choose A = 0.1. Clearly, it is possible to 
study a rich variety of interplays between the two quasiparticle 
bands in terms of different symmetries for the d xz - and d yz - 
orbitals and with different gap magnitudes. Here, however, 
our main aim is to investigate how the conductance spectra are 
influenced by the non-trivial Fermi surface topology and dis- 
persion relations, and see how this compares with the cylindri- 
cal/spherical Fermi-surface and free-particle dispersion sce- 
nario employed in the usual BTK-paradigm. In particular, 
this is relevant to the interpretation of the point-contact spec- 
troscopy measurements of Refs i 1Q i n . There is, however, an 
important cavaet with regard to which conclusion one may 
draw with regard to the symmetry of the superconducting OP 
from the tunneling data of Refs. 101 1 . In these works, polycrys- 
talline samples were used, while the orbital/nodal structure of 
the OP can only be convincingly probed in single crystal spec- 
imens. This is because tunneling into polycrystalline samples 
may lead to intrinsic averaging effects which distort the con- 
tribution from anisotropic OPs. 

In Fig. [2] we plot the conductance for tunneling along the 
(lOO)-direction for several OP symmetries. As seen, the d xy - 
wave case stands out from the rest as it features a considerable 
ZBCP. Comparing with the experimental data of Ref. [TH we 
would conclude that a nodal d-wave OP is likely to be realized 
in LaOo.gFo.i-^FeAs. The results of Ref. 10 seem to be most 
consistent with either s-wave or extended s-wave pairing, as 
only one gap is seen in the spectra. For the s-wave and d xy - 
wave cases, the standard BTK approach appears to suffice in 
order to qualitatively say something about the OP symmetry. 
However, the results are quite different from the usual BTK- 
approach when considering the extended s-wave and d x 2_ y 2- 
wave symmetries. More specifically, we find satellite features 
at subgap energies, including sharp peaks. These features 
most likely pertain to the specific band-structure which we 



consider here (see Fig. [TJ, and are thus not possible to capture 
within the conventional BTK-treatment with the cylindrical 
Fermi surface approximation. In fact, the DOS in our min- 
imal two-band model is a highly non-monotonic function of 
energy and contains two van Hove singularities. 817 . Let us 
also consider the case where there is one fully gapped OP and 
one nodal OP to see what fingerprints this combination leaves 
in the conductance spectra. In Fig. [3j we plot the conductance 
for the case where Ak x is fully gapped, while Ak y has a nodal 
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FIG. 3: (Color online) (a) and (b): Plot of the conductance spectra for 
tunneling along the (100)-axis in an iron-pnicitide N|S junction for 
the case of one fully gapped OP and one nodal OP. (c) Evolution of 
the Fermi surface topology for [i = {1.24, 1.54, 1.84} in the direc- 
tion of the arrow, (d) Conductance spectra for the c x + c y symmetry 
with Z — 3 for different doping levels. 
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symmetry. For concreteness, we consider s-wave + d^-wave 
pairing and d x 2_ y 2-wscve + d xy -wave pairing in Fig. [2 a) and 
(b), respectively. As seen, the nodal OP gives rise to a ZBCP 
while there are several satellite features in addition to the large 
coherence peak at the gap edge. The plots are qualitatively 
similar regardless of whether the fully gapped OP is s-wave or 
d x 2_ y 2-wave, while the features in the conductance are quali- 
tatively more pronounced in the s-wave case due to the better 
gapping of the Fermi surface. Finally, we consider the evo- 
lution of the conductance spectra upon changing the doping 
level fi. The Fermi surface topology evolves with a change in 
ji as shown in Fig. [3c): the electron-pockets increase in size 
while the hole-pockets decrease in size upon increasing fi. To 
see how the subgap features obtained in Ref. [2] evolve upon 
modifying fi, consider Fig. [3d) where we consider the c x +c y 
symmetry with Z = 3. As seen, the satellite features shown 
in Fig. [2] are still present and qualitatively the same, but they 
are shifted to different bias voltages. 

Summary. In summary, we have developed a theory for 
Andreev reflection in the iron-based high-T c superconduc- 
tors. Starting with a tight-binding model on a square lattice 



to model the puckered FeAs planes, we have investigated sev- 
eral order parameter (OP) symmetries and the resulting con- 
ductance spectra. Taking fully into account the Fermi surface 
topology and the quasiparticle dispersion relation, we have 
investigated scenarios where the symmetry of the supercon- 
ducting OP in both bands is the same and where it is differ- 
ent, i.e. one is fully gapped and the other is nodal. We find 
that the standard Blonder- Tinkham-Klapwijk (BTK) formal- 
ism should give qualitatively correct results for the case where 
the OP symmetries on both bands are either isotropic s-wave 
or ti-wave. However, the results differ considerably for the 
extended s-wave symmetries, as we find satellite features at 
subgap energies which are absent within the usual BTK treat- 
ment. Our results may be useful in the context of analyzing 
quantum transport data of tunneling in normal | superconductor 
junctions involving the iron-pnici tides. 
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